arXiv:1508.00678vl [astro-ph.EP] 4 Aug 2015 


Research in Astron. Astrophys. 2015 Vol. X No. XX, 000-000 
http://www. raa-journal. org http://www. iop. org/joumals/raa 


Research in 
Astronomy and 
Astrophysics 


Global Mg/Si and Al/Si Distributions on Lunar Surface Derived 
from Chang’E-2 X-ray Spectrometer * 

Wu-Dong Dong, Xiaoping Zhang, Meng-Hua Zhu, Aoao Xu, Zesheng Tang 

Lunar and Planetary Science Laboratory, Macau University of Science and Technology, Macao, China; 
Corresponding authors: xpzhang@must.edu.mo (X. Zhang) and dowudong@gmail.com (W.-D. Dong) 


Abstract X-ray fluorescence remote sensing technique plays a significant role in the 
chemical compositions research of the Moon. Here we describe the data analysis method 
for China’s Chang’E-2 X-ray spectrometer (CE2XRS) in detail and present the prelimi¬ 
nary results: the first global Mg/Si and Al/Si maps on the lunar surface. Our results show 
that the distributions of Mg/Si and Al/Si correlate well with the terrains of the Moon. The 
higher Mg/Si ratio corresponding to the mare regions while the lower value corresponding 
to the highland terrains. The map of Al/Si ratio shows a reverse relationship with the map 
of Mg/Si ratio. 
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1 INTRODUCTION 

Chemical compositions on the lunar surface are critical in the determination of geochemical nature of 
lunar terrains and the geologic evolution history of the Moon. X-ray spectroscopy is considered as one of 
the most effective method among all remote sensing techniques to study the major elements abundance 
on the airless planetary surface, such as the Moon, Mercury, and asteroids (Adler et al. 1973; Clark & 
Trombka 1997; Trombka et al. 2000; Grande et al. 2003; Sun et al. 2008; Nittler et al. 2011). In this 
technique, the solar X-rays interact with the planetary’s surface materials producing characteristic X-ray 
fluorescences at the uppmost surface layer. The characteristic fluorescent X-rays of major elements can 
be detected by the instrument in orbit, such as Ka lines of Mg (1.254 keV), Al (1.487 keV), Si (1.740 
keV), Ca (3.692 keV), Ti (4.511 keV), and Fe (6.404 keV). 

X-ray spectrometer was carried on many spacecrafts in the past. The X-ray spectrometer onboard 
the Russian Luna 12 orbiter sent in 1966 first successfully observed fluorescent X-rays from the Moon 
(Mandel’Shtam et al. 1968). Later in 1971 - 1972, X-ray fluorescences were measured by proportional 
counters in lunar equatorial region of the nearside in Apollo 15 and 16 missions, covering ~ 10 % of the 
global lunar surface. This observation was considered as the only observation of the large-scale chemi¬ 
cal compositions of lunar surface from the X-ray spectrometer. The Al/Si and Mg/Si ratios were derived 
from Apollo X-ray Spectrometer data, providing a preliminary understanding of the local distribution 
feature of the lunar surface chemical composition (Adler et al. 1972a,b; Clark 1979). The Demonstration 
of a Compact Imaging X-ray Spectrometer (D-CIXS) payload onboard SMART-1 launched by ESA 
in 2003 detected the characteristic X-ray lines from rock-forming elements during several solar flare 
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events (Grande et al. 2003; Grande et al. 2007; Swinyard et al. 2009). Unfortunately, the detector suf¬ 
fered severe radiation in orbit during its observation, from which no accurate quantitative analysis was 
presented. The X-ray spectrometer onboard Kaguya spacecraft that was launched in 2007 also suffered 
the same problem (Yamamoto et al. 2008; Okada et al. 2009). The X-ray spectrometers onboard the first 
Chinese lunar spacecraft Chang’E-1 launched in 2007 (Sun et al. 2008; Ouyang et al. 2008; Ouyang 
et al. 2010a,b) and Indian Chandrayaan-1 launched in 2008 (Grande et al. 2009) successfully detected 
the fluorescent X-rays from rock-forming elements (Ouyang et al. 2008; Peng 2009; Narendranath et al. 
2011; Weider et al. 2012, 2014). Since the Sun was in quiescent period at that time, the solar flares were 
few and the incident X-ray intensity was not powerful enough. As a result, the data is not sufficient for 
producing a global elemental distribution map. 

The Chang’E-2 (CE-2) spacecraft, which was the second unmanned China’s lunar probe, was suc¬ 
cessfully launched on 1 October 2010. An innovative X-ray spectrometer (CE2XRS) was onboard CE- 
2. Ban et al. (2014) derived the elemental abundances of Mg, Al, Si, Ca and Fe in the lunar Oceanus 
Procellarum region using CE2XRS data acquired during an M-class solar flare event. However, there 
is no global elements abundance map of the Moon from X-ray spectrometer until now. Fortunately, the 
Sun was in its active period in which a number of solar flare events were observed during the CE-2 
mission . After more than half year’s observation, large amounts of scientific data were obtained by 
CE2XRS. It is, therefore, the best opportunity to derive the first global chemical composition maps of 
major elements with CE2XRS. 

In this paper, we describe the data analysis procedure and present the global Mg/Si and Al/Si distri¬ 
butions on lunar surface derived from CE2XRS data. This is the first global map derived by the remote 
X-ray fluorescence spectroscopy on the Moon. Lunar global geochemical features are discussed based 
on the element distribution maps. Three major geological units on the Moon (lunar mare, highland and 
South Pole-Aitken basin) are identified. 

2 CHANG’E-2 X-RAY SPECTROMETER 

The X-ray spectrometer, one of the scientific payloads onboard CE-2 spacecraft, was aimed to obtain 
major rock-forming elements distributions (i.e.. Mg, Al, Si, Ca, Ti, and Fe) on the Moon. Since the 
detailed structure and specification of CE2XRS has been presented by Peng et al. (2009b) and Ban et al. 
(2014), we here only give a brief description. 

The CE2XRS was designed as a compact detector. It consisted of lunar X-ray detector, solar X-ray 
monitor and electric box (see Figure 1), which was an improved one compared to the X-ray spectrometer 
on Chang’E-1 spacecraft. A pair of orthogonal collimators made of permanent magnet limits the field of 
view of the detectors to make the best spatial resolution of ~ 200 km x 70 km at a distance of 100 km 
from lunar surface. The layout of the detectors and collimators were described by Peng et al. (2009a). 
The lunar X-ray detector was designed to measure the lunar fluorescent X-rays. It was composed of a 
soft X-ray detector array SXDA and a hard X-ray detector array (HXDA). The SXDA had 4 Si-PIN 
detectors named 111 1-9, 1111-10, 2222-9, 2222-10 . The detective energy range of SXDA was 0.5-10 
keV and the energy resolution was 300 eV@5.95 keV. As for HXDA, the detective energy range was 25- 
60 keV and the energy resolution was 6.0 keV@59.5 keV. The solar X-ray monitor (SXM) had a small 
Si-PIN detector . It was the same as the detector in SXDA except the effective detector area. The SXM 
was designed to measure the solar X-ray spectrum simultaneously. A 55 Fe radioactive source emplaced 
on the 1111-9 detector in SXDA was used to calibrate the instrument in orbit. X-rays emitted from the 
55 Fe source are the Mn Ko and K/3 lines at 5.898 keV and 6.490 keV, respectively. The instrument 
specifications of the X-ray spectrometer are listed in Table 1 . 

3 CE2XRS DATA ANALYSIS 

The CE-2 spacecraft kept a polar orbit with a period of ~ 118 minutes at an average altitude of 100 km. 
There were 2739 orbits data collected from October 15, 2010 to May 20, 2011. In our analysis, we use 
the CE2XRS level 2C dataset and only analyze the data from the detector 1111-9 in SXDA. The data 
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Fig-1 The CE2XRS instruments, from left to right listed are electronic box, solar X-ray 
monitor, lunar X-ray detector. The photograph is from Peng et al. (2009b). 


Table 1 Technical specification of CE2XRS from Peng et al. 
(2009b). 


Component 

SXD 

HXD 

SXM 

Detector 

Si-PIN (4 chips) 

Si-PIN (16 chips) 

Si-PIN (1 chip) 

Filter 

12.5 gm Be 

1 fim Al 

12.5 /im Be 

Effective area 

1 cm 2 

16 cm 2 

0.2 mm 2 

Detective range 

0.5-10 keV 

25-60 keV 

0.5-10 keV 

Energy resolution 

300 eV@5.9 keV 

6 keV@59.5 keV 

300 eV@5.9 keV 


format is shown in Table 2. The description of the dataset was presented by Ban et al. (2014). A general 
X-ray spectrometer data analysis method was published by Clark & Trombka (1997) , Starr et al. (2000) 
and Peng (2009). Here we describe the CE2XRS data analysis method for global elemental distribution 
map in detail. 


Table 2 The format of CE2XRS level 2C dataset. 


Column 

Label 

Column 

Label 

1 

time 

7 

instrument_azimuth_angle_2 

2 

longitude 

8 

solar_incidence_angle 

3 

latitude 

9 

detector_number 

4 

distance 

10 

energy 

5 

instrument_incidence_angle 

11 

counts 

6 

instrument.azimuth.angle 

12 

quality .state 


3.1 Data Check 

For the CE2XRS level 2C data, there is a label “quality.state” (see Table 2) appended to the end of each 
record. The value of “quality.state” equaled “00” (decimal value = 0) indicating that the record is good. 
If a record label “quality.state” is not equal to “00”, this record is eliminated (see Figure 2). The label 
“distance” is the altitude between the CE-2 spacecraft and the lunar surface. The normal average work 
altitude of CE2XRS is ~ 100 km. The records with the value of label “distance” greater than 125 km are 
deleted (see Figure 3), although the percentage of such measurements is very low. 
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Fig. 2 The scatter plot of label “quality_state” (QS) versus time from the data in orbit number 
578. The time means the number of days since October 1, 2010. This orbital record started 
from 11:22 to 13:20 on November 24, 2010 (UTC). The bad record with the “quality_state” 
“40” happened at 11:38 and 12:41. These abnormal records are eliminated in the data check 
process. 


Also, we check the relationship between channel number and energy for the detector 1111-9. 
This detector has 1024 channels. The channel information is stored in the level 2B dataset. The level 
2C dataset is generated when the channel-energy calibration was carried out on the level 2B dataset. 
This preliminary calibration was made by Institute of High Energy Physics and National Astronomical 
Observatories. We find a good linear relationship between channel number and energy (see Figure 4) 
which demonstrates the good quality of the CE2XRS data. After the data check is finished, the final 
corrected data is ready for spectral analysis. 

3.2 The Solar Activity and Data Selection 

The solar X-ray is the primary excitation source for the fluorescent X-rays generation from lunar surface. 
When the Sun is in active period, the solar X-ray flux increase significantly. The flux of higher energy 
X-rays may increase several orders of magnitudes and makes the solar X-ray spectra much harder during 
solar flare. Hence, fluorescent X-rays from lunar surface increases greatly and the signal-to-noise ratio 
also goes up in this period. When the Sun is in quiescent period, the fluorescent X-rays is much lower and 
much more data accumulation time is needed to achieve significant results. The relationship between the 
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time 


Fig. 3 The scatter plot of label “distance” versus time from the data in orbit number 169. 
The time means the number of days since October 1, 2010. This orbital record started from 
00:25 to 02:23 on October 22, 2010 (UTC). The “distance” (orbit altitude) sharply increased 
at 02:07. These records are removed while the CE-2 spacecraft might be doing orbital transfer 
at that time. 


variation in solar activity and the intensity of fluorescent X-rays had been derived by Clark & Trombka 
(1997). 

In our analysis, the solar plasma temperature ( T so i ar ) parameter is used to represent the intensity of 
solar activity. We adopt the solar plasma isothermal model to derive the solar plasma temperature. The 
solar X-ray flux data are taken from the Geostationary Operational Environmental Satellite (GOES). 
GOES measures the solar flux in two channels (1-8 A and 0.5-4 A, corresponding energy is 1.55- 
12.4 keV and 3.1-24.8 keV) simultaneously in every 3s. Thomas et al. (1985) and Garcia (1994) have 
published methods for calculating solar plasma temperature from GOES data. This model was used 
by Trombka et al. (2000) and Nittler et al. (2001) in the NEAR-Shoemaker XRS papers to interpret 
the elemental composition of asteroid 433 Eros and also adopted in CE-1 XRS data analysis by Peng 
(2009) . The relationship between solar plasma temperature and the solar flux ratio (SFR) was obtained 
by fitting method (Peng 2009). The SFR is defined by the ratio of flux in channel 1-8 A and flux in 
channel 0.5-4 A. When the SFR increases, the solar plasma temperature goes down. 

From GOES solar flare event report , we find that there were many solar flare events during the 
CE2XRS observation period from October 15, 2010 to May 20, 2011. For the solar flare event with 

1 http://www.ngdc.noaa.gov/stp/space-weather/solar-data/solar-features/solar-flares/x-rays/goes/ 
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Channel Number 


Fig. 4 The linear relationship between channel number and energy from detector 1111-9 data. 
The points from left to right are Mg(Ka), Al(Ko;), Si(Ka), Ca(Ka), Ca(K/3), MniKo) and Mn 
(K/3). 


level B and above, the footprint of CE2XRS on the Moon covered all of the lunar surface. The solar 
condition was good for mapping out major elements, such as Mg, Al, Si, on the entire lunar surface. 
To obtain statistically significant results for globally mapping, the X-ray spectra with the same solar 
condition and at the same location are accumulated. To increase the signal-to-noise ratio, only the data 
acquired in the sunlit side of the Moon with solar incidence angle less than 88° are accumulated. The 
final data are categorized by different solar plasma temperatures. We find the amount of data with the 
solar plasma temperature of 4 MK is sufficient for deriving global distribution map of Mg, Al and Si 
on the entire lunar surface. Therefore, this solar condition is selected for our analysis. Hereafter, all 
operations are done on the data at this solar condition. 

The lunar surface is partitioned into a series of equal area grids for mapping CE2XRS data. 
According to the space resolution of CE2XRS and the statistics of the data, the grid size is set to be 
20° x 20°. The partitioned grids start from the lunar equator to the polars with the same latitude inter¬ 
vals but different longitude intervals. The total number of grids is 114. All corrected data in the same 
grid and in the same solar condition are accumulated to form a spectra for further analysis. 
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3.3 Background Determination 

The X-rays detected by CE2XRS in orbit include not only the fluorescent X-rays from the lunar surface, 
but also the backgrounds which come from internal electronic noise of the instrument (Peng 2009), 
external cosmic ray induced background and scattered solar X-rays (Nittler et al. 2001). Ban et al. 
(2014) made a comparison of the electronic noise background at different temperatures. It was found 
that the electronic noise of the detector was fairly constant and change slightly with the detector’s 
temperature. The overall shape of the cosmic ray induced background spectra also tends to be constant 
over the timescale of a few hours, but the magnitude varies (Lim & Nittler 2009; Nittler et al. 2011; 
Narendranath et al. 2011). In our analysis, the internal electric noise and cosmic ray induced background 
are determined by summing over all measurements in the same solar condition and in the dark side of 
the Moon. This method for background determination was widely adopted by many investigators (Clark 
& Trombka 1997; Peng 2009; Nittler et al. 201 1; Narendranath et al. 2011; Weider et al. 2012; Ban et al. 
2014). Ban et al. (2014) used the dark side with solar incident angle larger than 90° as the background. 
To avoid illuminations of Sun light and zodiacal light into the detector directly near the terminator 
region, only the spectra with a solar X-ray incident angle greater than 120° is added for background 
determination in our analysis. In our case, the background variation with incident angle larger than 120° 
is more stable as a function of incident angle. And then, the background is normalized according to 
the accumulation time and subtracted directely in the spectra acquired in the sunlit side of the Moon. 
Figure 5 shows an example spectra observed by CE2XRS in sunlit side and in dark side with solar 
incident angle larger than 120° when T so i ar = 4 MK. The scattered solar X-ray background in the 
sunlit side has been calculated theoretically by Clark & Trombka (1997). Nittler et al. (2001) have used 
this method to predict scattered solar X-rays spectra at different solar plasma temperatures. The shape 
in low energy part of scattered solar X-ray spectra is approximately a Gaussian profile at low solar 
plasma temperature. Therefore we use a Gaussian function to model the spectra of scattered solar X- 
rays background. As shown in the following section, this function describes the low energy tail of the 
spectra well. 

3.4 Spectral Deconvolution and Elemental Abundance Mapping 

After the electronic noise and cosmic induced background subtracted, the residual spectra were decon¬ 
volved to obtain the flux of each major element, such as Mg, Al and Si. For most spectra, the fitting 
function is the sum of three normalized Gaussian functions for Mg, Al, Si and a Gaussian function for 
scattered solar X-ray background (see Figure 6). Each element intensity is determined by the best fit 
of the spectra with a minimum x' 2 method that takes into account the error bar of each data point. The 
fitting procedure is coded and implemented in ROOT software. One can also obtain the uncertainties 
of the X-ray flux of each element from the fitting results. Based on the fitting results, we derived X-rays 
intensity ratio of each element relative to Si. Si abundances only vary slightly on the lunar surface. In 
the analysis, ratioing to Si helps to remove matrix effects that come from compositional variations in the 
lunar regolith and from geometric corrections (Clark & Trombka 1997). Then, the elemental flux ratios 
are used to create global distribution maps on the lunar surface. 

4 RESULTS AND DISCUSSIONS 

After above data process steps are carried out, we obtain the first global Mg/Si and Al/Si maps on the 
Moon derived from CE2XRS level-2C data on the condition that the solar temperature is 4 MK. The 
result is normalized by global Mg/Si and Al/Si data derived from Gamma Ray Spectrometer on Lunar 
Prospector (LPGRS) (Prettyman et al. 2006). Figure 7 shows the global Mg/Si map with a resolution 
of 20° x 20° and the same resolution map from LPGRS. Figure 8 shows the global Al/Si map with a 
resolution of 20° x 20° and the same resolution map from LPGRS. From Figure 7, we can find that the 
high Mg/Si ratios are mainly concentrated at the PKT (Procellarum KREEP Terrane) and SPA (South 
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Fig. 5 An example spectra observed by CE2XRS in sunlit side (in black point) and in the 
dark side with solar incident angle larger than 120° (in red point) when T so i ar = 4 M K. 


Pole-Aitken) while FHT (Feldspathic Highlands Terrane) has a low value. However, the high Al/Si ratios 
are inversely concentrated at the FHT and with a low value at PKT and SPA, as shown in Figure 8. All 
these maps are consistent well with the maps derived from LPGRS. 

The possible reasons of the differences in certain regions of elemental abundance ratios derived 
between CE2XRS and LPGRS are presented as follows. The solar spectra used for CE2XRS derivation 
in initial analysis is from GOES data, which may result in uncertainty to some extent. The solar activity 
calibration model that use CE2XRS Solar X-ray Monitor data will be established in future work. On the 
other hand, the detection depths of the two techniques are different. X-ray fluorescence is exicited by 
incident solar X-rays at the upmost surface of the lunar regolith. The penetration depth of solar X-rays 
is ~ 100 gm, while the penetration depth of gamma rays is on centimeter scale. 

5 CONCLUSIONS 

We have described data analysis procedures in detail for CE2XRS data and have presented the global 
distributions of Mg/Si and Al/Si on the lunar surface from CE2XRS, which are the first global maps 
derived from X-ray spectroscopy. Three major geological units (lunar mare, highland and South Pole- 
Aitken basin) are identified. The results demonstrate that X-ray remote-sensing technique can also pro¬ 
vide important compositional information about lunar surface. 
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Fig. 6 The spectra with electronic noise and cosmic induced background substracted are fitted 
by a sum of Gaussian functions for Ka lines of Mg (green line), A1 (blue line), and Si (pink 
line), and superimposed on a scatter solar X-ray background in yellow. The line in red is the 
best-fit result. 
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Fig. 8 The equal-area (20° x 20°) global map of Al/Si molar ratios distribution on lunar 
surface (Up: derived from CE2XRS; Down: derived from LPGRS). 
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